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We study wave propagation in strongly nonlinear ID diatomic granular crystals under an impact 
load. Depending on the mass ratio of the ‘light’ to ‘heavy’ beads, this system exhibits rich wave dy¬ 
namics from highly localized traveling waves to highly dispersive waves featuring strong attenuation. 
We experimentally demonstrate the nonlinear resonant and anti-resonant interactions of particles 
and verify that the nonlinear resonance results in strong wave attenuation, leading to highly efficient 
nonlinear energy cascading without relying on material damping. In this process, mechanical energy 
is transferred from low to high frequencies, while propagating waves emerge in both ordered and 
chaotic waveforms via a distinctive spatial cascading. This energy transfer mechanism from lower 
to higher frequencies and wavenumbers is of particular significance towards the design of novel 
nonlinear acoustic metamaterials with inherently passive energy redistribution properties. 

PACS numbers: 45.70.-n 05.45.-a 46.40.Cd 


Introduction. Granular crystals, i.e., periodically 
packed arrays of solid particles interacting elastically via 
nonlinear contacts, have recently become popular test 
beds for understanding fundamental structures emerging 
within nonlinear wave dynamics, such as traveling waves 
[1-3], shock waves [4-6], discrete breathers [7, 8], and 
nanoptera [9]. In these systems, tunable energy transport 
is of particular interest. For example, the wave propaga¬ 
tion speed can vary significantly from linear to highly 
nonlinear regimes by changing the external compression 
applied to the granular system [1]. The energy propaga¬ 
tion patterns can be tuned from highly localized to mod¬ 
ulated or widely dispersed shapes by introducing local 
resonances in constituents [9]. The controllable degree of 
disorder also changes the nature of energy attenuation, 
enabling both exponential and power-law type scenar¬ 
ios [10, 11]. This wealth of energy transport character¬ 
istics stems from the widely tunable interplay of nonlin¬ 
earity, discreteness, and heterogeneity in such media. 

Among the different heterogeneous variants of granu¬ 
lar crystals, an ordered diatomic granular crystal - in the 
form of alternating heavy and light particles - is one of 
the most prototypical ones [3, 12, 13]. Remarkably, this 
simple dimer lattice shows highly intriguing energy trans¬ 
port phenomena. Previous studies demonstrated that 
diatomic granular chains support various types of highly 
localized traveling waves, nonlinear beating pulses, and 
highly dispersive waves, depending on the mass ratio of 
the heavy and light beads [3, 12-14]. It is of particu¬ 
lar interest that wave attenuation can be maximized at 
specific mass ratios, while efficient transmission of en¬ 
ergy without attenuation can also occur at another set of 
mass ratios, resulting in highly nonlinear solitary travel¬ 
ing waves. This pertains to the so-called nonlinear reso¬ 
nance and anti-resonance mechanism [3, 13]. The mass 


ratios satisfying the resonance and anti-resonance condi¬ 
tions appear in an interlaced manner as the relevant pa¬ 
rameter is varied. Despite previous theoretical/numerical 
studies on these resonance and anti-resonance mecha¬ 
nisms, corresponding experimental work has been ex¬ 
tremely limited [15]. 

Here, we push considerably forth the experimental 
state of the art of the system by enabling a full-field 
visualization of lattice’s dynamics via laser Doppler vi- 
brometry. These experimental advances, in turn, lead us 
to consider crucial theoretical/numerical aspects, such as 
the nonlinear energy transfer mechanisms across length 
and frequency scales. Specifically, we investigate how the 
resonances chiefly disintegrate and transfer the energy 
into stable periodic modes, leading to extremely efficient 
nonlinear energy scattering and impact mitigation mech¬ 
anisms. The remainder of the energy is partitioned to 
a wide range of frequencies and modes, resulting in an 
apparently chaotic spatial tail dynamics. This set of fea¬ 
tures is in sharp contradistinction with the anti-resonance 
scenario where the energy appears to localize in traveling 
wave quanta (i.e., isolated wavepackets). 

Our approach towards exploring the energy transfer 
mechanisms in granular crystals, in addition to substan¬ 
tial experimental developments, is motivated theoreti¬ 
cally by the usefulness of this type of approaches in re¬ 
vealing cascading dynamics in a large variety of other sys¬ 
tems. These include classical and quantum fluids [17, 18], 
plasmas [19], and solid structures such as vibrating thin 
plates [20-22]. Previous studies tackled the topic of en¬ 
ergy distribution in granular crystals in the context of ta¬ 
pered and decorated chains [2, 5], disordered settings [4], 
chains bearing impurities [23], interaction of waves with 
their own radiation [24], and that of narrow or broad¬ 
band frequency signals with waves stemming from reflec- 
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FIG. 1: Schematic of experimental setup for measuring non¬ 
linear wave propagation in a diatomic granular crystal. Inset 
shows an instrumented sensor particle. 


tions [25]. Nevertheless, the mechanism of energy cascad¬ 
ing across time- and length-scales has not been suitably 
explored yet. We believe that this study will plant the 
seeds towards the investigations of such energy cascading 
dynamics in granular crystals as well as nonlinear lattice 
systems more generally. 

Experimental and numerical setup. Figure 1 shows the 
test setup consisting of a granular chain, a striker, and a 
laser Doppler vibrometer (LDV). The chain is composed 
of alternating ‘heavy’ and ‘light’ spherical particles. All 
particles are made of chrome steel (E = 210 GPa, v = 
0.29, p = 7,810 kg/m 3 ). The radius of the heavy particles 
is R = 9.525 mm. We alter the radius of the light particles 
to set various mass ratios for the chain (see Supplemen¬ 
tal Document [24] for details). We position the smaller 
particles on polytetrafluoroethylene (PTFE) rings tied to 
stainless steel rods to align the center of mass of all par¬ 
ticles (Fig. 1). We generate nonlinear stress waves in the 
granular crystal by impacting one end of the chain with 
a striker, which is identical to the heavy particle. The 
striker is released from a ramp and attains an impact ve¬ 
locity of 0.94 dz 0.017 m/s, hitting the first bead which 
is a heavy bead with an embedded piezoelectric ceramic 
disk at its center (see the inset of Fig. 1). At the time 
of impact, this sensor bead generates high voltage, which 
triggers our data acquisition system. The other end of 
the chain is fixed with a heavy steel wall. 

We position the LDV at a slanted angle (45° in this 
study) to measure the velocity of a particle in each im¬ 
pact event (Fig. 1). We repeatedly measure the parti¬ 
cle’s velocity by shifting the LDV along the guiding rail 
and synchronize the measured data with respect to the 
striker’s impact moment to reconstruct the wave prop¬ 
agation profiles. We examine the propagation of waves 
in two anti-resonance cases, labeled AR1 and AR2 with 
mass ratios of 1 (the homogeneous case, operating as a 
benchmark) and 0.34 respectively, and for two resonance 
cases, labeled R1 and R2 with mass ratios of 0.59 and 
0.24 respectively [3, 13]. 

We complement these experimental results by numeri¬ 


cal simulations with a discrete element model. The chain 
is modeled with point masses, which are connected by 
nonlinear springs representing the Hertzian contact be¬ 
tween spheres [26]. Therefore, the equations of motion of 
the diatomic granular chain can be written as follows: 

o/o q/ o 

M2n-lU2n-l = P[^2n-2 ~ “ P[u2n-l ~ ^2n]+ , 

M 2n U2n = P[u2n-1 ~ «2n]+ 2 “ /?[«2n “ «2n+l]+ 2 - 

Here M and u represent the mass and displacement of 
particles (odd and even subscripts correspond to heavy 
and light beads), (3 is the contact coefficient between two 
particles [27], and [s] + =max(s, 0), implying that the sys¬ 
tem does not support tensile force. We note that the 
dynamics of this system is fully re-scalable with energy, 
so for sufficiently small impulsive excitations the results 
do not change qualitatively [3, 13]. 

Results and discussion. Figure 2 shows numerical and 
experimental surface maps of the wave propagation in 
AR2, Rl, and R2. In Fig. 2(a), the impact leads to the 
formation of a number of robust solitary pulses (three 
of these pulses are directly visible), which quantize the 
energy repartitioning it primarily into the first one [27- 
29]. The experimental results in Fig. 2(d) obtained by 
the LDV corroborate the numerical simulations, though 
the tertiary packet is not distinguishable due to the dis¬ 
sipation in experiments. It should also be noted that 
Figs. 2(a) and (d) employ different time and length scales 
because of the shorter granular chain used in experi¬ 
ments. In this anti-resonance case, the traveling waves 
propagate with almost constant speeds without shedding 
energy behind the leading pulses. The velocities of these 
wave packets follow the relationship, V ~ i^/ 6 , where 
F m is the peak amplitude of the dynamic force. This is 
the same as the conventional solitary wave in a homoge¬ 
neous chain (i.e., AR1), which shows a single hump with 
its wavelength of ^5 particles’ diameters [1]. In AR2, the 
wavelength is wider than that in AR1 (see Supplemental 
Material [27]), and the heavy and light particles exhibit 
different velocity profiles [3]. We experimentally verified 
the waveform in AR2 as shown in the inset of Fig. 2(d). 

In sharp contrast to anti-resonances, resonances cause 
the propagating wave to shed its energy behind in the 
form of oscillating tails, thereby attenuating the leading 
pulses (see Figs. 2(b) and (c) for Rl and R2, respectively) 
[13]. We accurately measured the particle velocity pro¬ 
files in these cases as well (see Figs. 2(e) and (f)). In 
Rl, the primary wave loses a considerable fraction of its 
energy by forming a wave tail oscillating initially in a 
regular fashion. However, after approximately 80 par¬ 
ticles, three leading peaks form a wave packet (see the 
inset in Fig. 2(b)). These pulses exchange energy among 
them as the wave packet propagates in the medium. This 
behavior is attributed to the so-called nonlinear beating 
phenomenon [13]. We observe that once the nonlinear 
beating starts, the energy transfer from the leading wave 
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FIG. 2: Numerical and experimental surface maps for (a, d) AR2, (b, e) Rl, and (c, f) R2. The insets in (b) and (c) show 
an enlarged view of nonlinear wave beating phenomenon in Rl, and regular (I) and irregular (II) forms of wave tail in R2, 
respectively. The insets in (d-f) illustrate comparisons of three successive (heavy-light-heavy) particles’ velocity profiles between 
numerical (solid blue(heavy) and green(light)) and experimental (dashed red) results. 


packet to the wave tail is substantially reduced, as the 
relevant waveform detaches itself from the rest of the 
radiative tail. In R2, however, the primary wave con¬ 
tinuously loses its energy to the wave tail. In this case, 
as well as in that of higher order resonances, the oscilla¬ 
tions of the wave tail retain a very regular structure just 
behind the primary wave. They eventually transform to 
apparently chaotic forms within the far-held [see region 
(I) and (II) in the inset in Fig. 2(c)]. 

To further investigate the energy dispersion mecha¬ 
nism, we consider the relevant phenomenology in the 
Fourier domain. Figure 3(a) shows the frequency spec¬ 
trum of the particles’ velocity profiles in AR2, where the 
majority of energy is concentrated in low frequencies. In 
Rl (Fig. 3(b)), however, we find that a strong high fre¬ 
quency signal (~9 kHz) appears, and the magnitude of 
the low frequency signal decreases during the primary 
pulse’s shedding its energy to the wave tail (around up 
to 80 beads). This predominant mechanism of energy 
transfer ceases to exist as soon as the nonlinear beating 
starts. This represents the fact that the mechanical en¬ 
ergy carried initially by the low-frequency primary pulse 
is redistributed to the higher-frequency contained in the 
wave tail due to the nonlinear resonance mechanism, but 
is subsequently trapped in the beating mechanism with¬ 
out further substantial decay. In R2, this energy transfer 
from low to high frequencies happens persistently (see 
Fig. 3(c)) owing to the continuous energy shedding to the 
wave tail. We can discern in the relevant frequency pat¬ 
terns the presence of regular oscillations in the wave tail 
[corresponding to part (I) in the inset of Fig. 2(c)]. Ad¬ 


ditionally, there exists a fraction of the energy deposited 
to a wide range of frequencies corresponding to a dis¬ 
persion of the energy in other modes, resulting in the 
apparently chaotic oscillations of section (II) within the 
inset of Fig. 2(c). 

Equally telling as the frequency-domain responses, are 
the wavenumber plots in Figs. 3(d-f). In contrast to AR2, 
the resonance cases reveal how the energy of the primary 
wave at small wavenumbers is converted to a wavenumber 
around 7r/2 in the wave tail [see Figs. 3(e) and (f) in com¬ 
parison to Fig. 3(d)]. Remarkably, this length scale corre¬ 
sponds to the length of two sets of heavy and light beads. 
We note that the two strong signals near it/ 2 in R2 are 
because of the Fourier transform of mixed signals from 
the heavy and light beads in an alternating pattern (see 
Supplemental Document [27]). This wavenumber trend 
suggests the emergence of a definitive periodic traveling 
wave pattern that is being excited in association with 
a periodicity of two-beads. Such periodic patterns are 
witnessed experimentally in Figs. 3(e) and (f), and they 
have also been studied analytically in dimer chains [14]. 
Our wavenumber analysis reveals the excitation of these 
states in a transient way for Rl, until the beating pat¬ 
tern forms. In R2, contrarily, the pattern is persistent 
and clearly discernible in region (I) of Fig. 2(c), enabling 
the continuous transfer of energy from the principal wave 
to the associated background. Once again in R2, when 
the wave tail transforms to chaotic oscillations, the length 
scale also spreads out widely as denoted by the speckles 
in Fig. 3(f). 

An alternative diagnostic towards the characterization 
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Particle Number Time [ms] 

FIG. 3: Frequency and wavenumber spectra of particles’ ve¬ 
locity profiles in (a, d) AR2, (b, e) Rl, and (c, f) R2, obtained 
from the analysis of the computational results and illustrating 
the presence of cascades associated with the redistribution of 
energy among different scales. The wavenumber is calculated 
based on the wavelength normalized by the average diameter 
of the large and small particles. 


of the energy transmission is provided in Fig. 4 under var¬ 
ious anti-resonance and resonance conditions. Here, we 
compare the energy carried by the leading wave packet 
after normalizing it with respect to the impact energy in¬ 
jected by the striker [27]. We observe highly efficient and 
constant energy transmission in AR1 and AR2 due to 
the formation of solitary pulses. The normalized energy 
is less than 1 in AR2, which is caused by the trailing soli¬ 
tary waves generated in the transition region as explained 
above. In sharp contrast to such anti-resonance cases, the 
primary pulse experiences a fast and dramatic attenua¬ 
tion in Rl. About 80% of the energy is transferred to the 
wave tail within just 3 ms, which corresponds to the time 
the wave takes to propagate through the first 80 particles 
(see Fig. 2(b)). Hence, the emerging periodic traveling 
wave is extremely efficient in rapidly draining the energy 
of the leading pulse. However, once the nonlinear beat¬ 
ing pulse is formed, the transfer of energy to the wave 
tail is essentially suppressed. On the other hand, in R2, 
the reduction of normalized energy is slower compared to 
that in Rl. Nonetheless, such wave attenuation happens 
continuously given the persistence in the latter case of 
the periodic traveling wave. As a result, the normalized 
energy carried by the primary pulse approaches asymp¬ 
totically to zero. This eventually leads the attenuation 
of the primary pulse in R2 to become larger than that in 
Rl. Higher order resonances share the principal charac¬ 
teristics of R2, as concerns the above phenomenology. 


The distinctive energy scattering phenomena between 
Rl and R2 stem from the difference in relative motions 
of the heavy and light beads. When the particles are 
squeezed by the leading wave, light particles move faster 
than heavy particles. The disparity between the heavy 
and light particles’ motions becomes more evident as the 
mass ratio decreases. In Rl, the light and heavy beads 
move at relatively similar frequencies due to the small dif¬ 
ference in their masses (mass ratio of 0.59). This results 
in the fast energy scattering under the facilitated transi¬ 
tion of energy between the heavy and light particles. In 
R2, however, the light particles’ motions are much faster 
than those of heavy particles due to the small mass ratio 
(only 0.24). This results in a slow - but more consistent 

- reduction of energy, leading to the higher wave attenu¬ 
ation performance than Rl. Given the limited length of 
the granular crystal tested in this study, we obtained ex¬ 
perimental data up to ^0.8 ms (inset of Fig. 4). We find 
that the experimental results corroborate the numerical 
simulations. 

These observations provide us with a complete picture 

- in both real and Fourier space, in both time and fre¬ 
quency domains - of how wave attenuation mechanisms 
can be achieved on the basis of nonlinear resonance mech¬ 
anisms. This type of process is fundamentally different 
from the conventional energy attenuation mechanisms re¬ 
lying on material damping and/or structural deformation 
effects. Here, the wave attenuation is achieved by max¬ 
imizing the energy dispersion and redistribution within 
the chain of granular particles. The existence of damp¬ 
ing in this system only enhances the wave attenuation. 
This is confirmed by our experimental data, where the 
energy of the primary pulse attenuates more than that 
in numerical simulations (Fig. 4). Moreover, the energy 
cascading from low to high frequency and large to small 
length scale plays a crucial role in a way partially rem¬ 
iniscent of similar cascading phenomena within fluidic 
systems. 

Conclusions. In the present work, we experi¬ 
mentally and numerically investigated the nonlinear 
resonance/anti-resonance phenomena in ordered di¬ 
atomic granular crystals, and their connection to highly 
effective nonlinear frequency redistribution and cascad¬ 
ing of energy across different scales. We verified the ex¬ 
istence of efficient energy transfer mechanisms for res¬ 
onances, and robust traveling solitary waves for anti¬ 
resonances. In particular, we demonstrated that the non¬ 
linear resonance mechanism can be highly useful in at¬ 
tenuating impact energy by dispersing it to the wave tail. 
In this process, we found an interesting energy transition 
mechanism, i.e., the low-to-high frequency energy trans¬ 
fer and large-to-small scale energy cascading. These cas¬ 
cades resemble analogous cascades encountered in turbu¬ 
lent flows, which shows that acoustic metamaterials such 
as the granular diatomic chain considered herein can be 
tuned to possess this type of turbulence-like behavior. 
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FIG. 4: Normalized energy carried by the primary wave 
packet in anti-resonances and resonances. Inset shows the 
comparison between the experimental and numerical results. 
The error bars represent standard deviations based on 5 tests. 

A dominant role within this mechanism was revealed to 
be played by periodic traveling waves which appear to 
be excited either transiently (Rl) or permanently (R2). 
This can offer a new way to designing nonlinear acoustic 
metamaterials based on inherently passive energy redis¬ 
tribution principles for impact mitigation purposes. 
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SUPPLEMENTAL MATERIAL: 

Nonlinear Low-to-High Frequency Energy Cascades in Diatomic Granular Crystals 


PACS numbers: 45.70.-n 05.45.-a 46.40.Cd 


I. HERTZIAN CONTACT AND TOTAL ENERGY OF EACH PARTICLE 

In the diatomic granular chain, contact interaction between two beads can be expressed by the Hertzian law as 
follows: 


F = p[u n - u nU ] 3 f (1) 

where /3 is a contact coefficient which depends on material properties and geometry of the two objects in contact. In 
case of a contact between large and small beads made of an identical material, the coefficient is defined as follows: 


2 E I RlRs 

3(1 -v*)\ R l + Rs 


( 2 ) 


where E and v are Young’s modulus and Poisson’s ratio of the beads respectively, and Rl and Rs are radius of the 
large and small beads, respectively [1]. In the case of the last bead with radius R is in contact with flat surface of the 
wall made of identical materials with the bead, its contact coefficient becomes: 


2 EVR 
3(1 -v 2 ) 


(3) 


The equation of motions of the granular chain shown in the manuscript are solved numerically using the Runge- 
Kutta method with a 10 -7 s time step. For this numerical simulation, we use a 300-particle chain -although only a 
fraction of these, typically 150 are shown in simulations e.g. of Fig. 2-. The numerical simulation shows accuracy 
of about 0.004% error in terms of the total energy conserved during a 20 ms period. It should be noted that we do 
not consider damping effects in the numerical simulation in order to focus on the wave attenuation caused solely by 
energy re-distribution in the chain. If we assume that the potential energy between the two particles in contact is 
equally distributed to the two particles, the total energy of the n-th particle at time t can be defined as follows: 


E n (t) = 7 }M n ii 2 n + ^/3[u ra _ 1 - u n fl 2 + - u„+i]+ 2 (4) 

To calculate the energy of the leading wave in Fig. 4 of the main manuscript, we conduct the summation of the energy 
for seven particles, which are wide enough to cover the envelop of the leading pulse in AR1, AR2, and R2 cases. In 
Rl, however, we add the energy of nine particles to account for the total energy of the beating waveform. The total 
energy of the leading waves is normalized with respect to the initial excitation energy, Eq = ^MoUq, where Vo is the 
impact velocity of a striker. 


II. PARAMETERS OF GRANULAR CRYSTALS IN EXPERIMENT 

In experiment, we measure wave propagation in two anti-resonance models, AR1 and AR2, and two resonance 
models, Rl and R2. Their dimensions and mass ratios are summarized in Table I. 


III. WAVENUMBER OF VELOCITY PROFILES IN THE FIRST (Rl) AND SECOND (R2) 

RESONANCES 


The high wavenumber components in the nonlinear resonance Rl and R2, shown in Fig. 3(e) and 3(f) respectively, 
correspond to the regular wave tail discussed in Fig. 2 (b) and (c), where ‘heavy’ and ‘light’ beads oscillate in an 
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FIG. 5: Wavenumber analysis in the case that two different profiles with an identical wavenumber of 7t/2 are mixed together in 
an alternating way. Based on the FFT analysis, (b), (d), and (f) show wavenumber of the wave profiles in spatial domains for 
(a) heavy beads, (c) light beads, and (e) the mixed data of the two, respectively. In (a) and (c), red curves represent continuous 
signals, while blue circles denote the sampled data at particle positions. 
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FIG. 6: Wavenumber analysis in the case that two different profiles with an identical wavenumber of 0.487T are mixed together 
in an alternating way. Based on the FFT analysis, (b), (d), and (f) show wavenumber of the wave profiles in spatial domains for 
(a) heavy beads, (c) light beads, and (e) the mixed data of the two, respectively. In (a) and (c), red curves represent continuous 
signals, while blue circles denote the sampled data at particle positions. 


TABLE I: Mass ratio and chain length of each experiment models 


Model ID 

Large bead 
diameter [mm] 

Small bead 
diameter [mm] 

Mass ratio 

Chain length [mm] 

(number of beads) 

AR1 

19.050 

19.050 

1 

609.6 (32) 

R1 

19.050 

15.978 

0.59 

595.5 (34) 

AR2 

19.050 

13.333 

0.34 

615.3 (38) 

R2 

19.050 

11.839 

0.24 

605.9 (39) 
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ordered fashion in the near field of the primary pulse. When we perform Fast-Fourier-Transformation (FFT) of velocity 
profiles of all beads (i.e., both heavy and light ones) in the spatial domain, we observe one or more wavenumbers 
symmetric with respect to 7r/2, which corresponds to the wave length identical to the 4 particles’ distance. Here, 
in the calculation of wavenumber we use the number of particles, instead of their actual dimensions for the sake of 
simplicity. The symmetry with respect to tt/2 is due to the FFT of a mixed signal composed of two different velocity 
profiles. For example, the heavy particles show small amplitude oscillations, while the light particles show larger 
amplitude ones due to the smaller inertia. If we separate the signals into two profiles, ‘heavy’ particle’s velocity profile 
and ‘light’ particle’s velocity profile (see Fig. 5 (a) and (c)), the FFT of these signals show wavenumbers up to i r, 
because the lattice constant is two particles in this case (see Fig. 5 (b) and (d)). These signals are symmetric with 
respect to the central wavenumber 7r/2 (see Fig. 5 (b) and (d)) due to the nature of the FFT of real valued-signals. 
This is also referred to as a spatial aliasing effect. It is noted that one is mirrored image of the other. Therefore, if 
the wavenumber of the harmonic signal is 7r/2, we only can see a single peak at tt/2 (Fig. 5(b) and (d)). Otherwise 
the FFT shows two peaks in symmetry with respect to tt/2 (Fig. 6 (b) and (d)). 

In the original signals where the two velocity profiles are mixed in an alternating pattern, the lattice constant is a 
single particle. Therefore, the FFT of this signal shows wavenumbers up to 27r, and it is symmetric with respect to 7r. 
Interestingly, in this FFT signal of the combined data, the symmetric wavenumbers in the separated velocity profiles 
(in the range of 0 to i r) appear on the top of each other in the same wavenumber range, from 0 to 7r. Therefore, when 
the oscillations of both ‘heavy’ and ‘light’ particles have a single wavenumber of 7t/2, only one wavenumber (tt/2) 
appears in the FFT of the combined signal within the range from 0 to 7r. On the other hand, when they have a same 
wavenumber which is other than 7 t/ 2, double peaks symmetric with respect to 7 t/ 2 appear. If the two signals have 
different wavenumbers, the FFT signals show four peaks in symmetry with respect to tt/2. 

In the first resonance (Rl), the high wavenumber component is initially dispersed, and then it is focused around 
7r/2 as the wave propagates up to about 80 particles. In the second resonance R2, however, two distinctive lines 
continuously appear as the wave propagates, which represent that the wavenumbers of the heavy and light particles’ 
oscillations are similar to each other with a minor offset from the tt/2. Again, this difference between Rl and R2 
stems from the difference of relative particles’ motions as we explained in the main manuscript. 


IV. EFFECT OF STRIKER MASS ON THE FORMATION OF SOLITARY WAVE TRAINS 

In a homogenous chain (AR1), a solitary wave train appears when the striker mass is larger than the particle mass 
in the chain [3, 4]. However, in a diatomic chain having the mass ratio of light to heavy particle less than 1 (e.g., 
AR2), solitary wave train always appears regardless of the striker mass. This is because the striker impact results in 
multiple collisions among light and heavy particles at the beginning of the chain, which, in turn, generates multiple 
solitary waves. Figure 7 (a) shows peak forces of the solitary wave train in AR2 at various mass ratios of the striker 
mass (M 0 ) to the heavy mass of the chain (M). Here we only show the peak forces larger than 0.2 N. In this chain, the 
first particle is a heavy bead, and the mass ratio of the light bead (m) to the heavy bead (M) is 0.34. We observe that 
even when the striker mass is less than the light mass (Mo/M = 0.17), a solitary wave train is formed and propagates 
as shown in Fig. 7(b). The inset shows the solitary wave train formed in the transition region. From Fig. 7(a), 
we also notice that when Mq/M is larger than a unity, the number of peaks in the solitary wave train substantially 
increases. 


V. WAVELENGTH OF THE SOLITARY WAVES IN DIATOMIC GRANULAR CRYSTALS 

In this section, we compare the wavelength of the solitary waves in diatomic granular chains at various mass ratios. 
We obtain the wavelength values from three different approaches: (1) the analytical method by Porter et al. [2] (solid 
blue curve); (2) the numerical simulations via a discrete element method (green squares); and (3) the experimental 
results via laser Doppler vibrometry (red circules) as shown in Fig. 8. In the numerical simulation and experiments, we 
count the number of particles (n) having a velocity larger than a threshold value in a fixed time frame, and determine 
n + 1 as the wavelength considering the smoothly decreasing wave edges. Here, the threshold is set to be 0.1% of the 
maximum particle’s velocity. It is to be noted that the wavelength varies depending on the measurement moment 
due to the discreteness of the system, and it can be overestimated at a small threshold value. Therefore, we process 
multiple data points at various time frames for statistical analysis. In the analytic expression, the wavelength is 
about five particles’ diameters at mass ratio 1, and it gradually increases and approaches approximately ten particles’ 
diameters as the mass ratio approaches zero. In the experimental and numerical results, the wavelength in AR1 is 
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FIG. 7: Solitary wave train in AR2 at various striker masses (a) magnitude of the peak contact forces for different striker 
masses (Mo) (b) image map of a solitary wave train in AR2 when the striker mass is half of the light mass (i.e., Mq/M — 0.17). 



FIG. 8: Wavelength of the solitary waves in a diatomic chain at various mass ratios from AR1 to AR8. The solid curve 
represents analytical prediction based on Porter et al. [2]. Green squares and red circles are from numerical simulations and 
experiments, respectively. 


slightly larger than five particle’s diameters. It is generally known that the wavelength of the conventional solitary 
wave (in AR1) is about 5 particles’ diameters. However, it has also been reported that the wavelength in AR1 is close 
to 7 particles’ diameters, if we account for extremely small wave edges [5]. This can be captured when we lower the 
threshold value of detectable wave velocity significantly. 
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